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We investigate the phenomenon of Bose-Einstein condensation on manifolds constructed 
as a product of a three-dimensional Euclidian space and a general smooth, compact d- 
dimensional manifold possibly with boundary. By using spectral (-function methods, we are 
able to explicitly provide thermodynamical quantities like the critical temperature and the 
specific heat when the gas of bosons is confined, in the three-dimensional manifold, by the 
experimentally relevant anisotropic harmonic oscillator potential. 
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Under specific circumstances a gas of spin-0 particles (bosons) undergoes a process of phase 



transition where the particles tend to reside in the lowest energy state (see, e.g. 
the well-known phenomenon of Bose-Einstein condensation la, Il3j which has become increasing' 

u n u 

interesting in the last few years because of its experimental realization; see, e.g., [2|, |6|, [9|, 



331). This is 

y 



2_j]. 



For the description of thermodynamical properties in these experiments, the case in w 



lich the 



Bose gas is confined by a harmonic oscillator potential is the most relevant one 



0, 0, 0,0, 



24 



27]. However, many other configurations have also been considered. For example, Bose-Einstein 



condensation has been studied in flat Minkowski space for rectangular enclosures 20l. l22l. l36l. l38l] and 
for more general, arbitrarily shaped, cavities [___]]. Also, more general confining potentials than the 



harmonic oscillator potential have been analyzed in 



29| . Generalizations of these investigations 



to curved spaces have been considered. In particular studies of Bose-Einstein condensation have 



)een performed for static Einstein manifolds in 



ay, 



4fJ] and for higher-dimensional spheres in 



33 ] . Moreover, Bose-Einstein condensation as a symmetry breaking phenomenon has been studied 



on static curved spacetimes of arbitrary spatial sections with and without boundary in 
also 



41 



12 



3J. 



43J; see 



In this work we utilize (-function regularization techniques, related to the ones developed in 
111 ] , in order to analyze the phenomenon of Bose-Einstein condensation on very general manifolds 
constructed as a product of the 3-dimensional Minkowski space and a d-dimensional smooth, com- 



* Electronic address: Guglielmo_Fucci@Baylor.edu 
' Electronic address: Klaus_Kirsten@Baylor.edu 



2 



pact manifold with or without boundary. More specifically, in the 3-dimensional Minkowski space 
we assume a confining harmonic oscillator potential as it is used in experiments. In particular, the 
dependence of the critical temperature and the specific heat on the dimension and the geometry 
and topology of the additional compact manifold is explicitly obtained. 

The main motivation for these studies is to understand how significant the impact of additional 
dimensions in particular on the critical temperature is. Given that the critical temperature can 
be determined to high accuracy, Bose-Einstein condensation experiments could provide a window 
into the world of extra dimensions by comparing experimental data with theoretical predictions in 
the presence of extra dimensions. This procedure has been successfully applied for example in the 
context of the Casimir effect and information about number and size of extra dimensions could be 

us 



established; see, e.g., 



13, 



251 ] . It is the aim of our article to start research in this direction. 
The outline of this article is as follows. In Section II we use heat kernel and zeta function tech- 
niques to analyze the partition sum of the system considered. In particular, the high temperature 
expansion of the partition sum in terms of heat kernel coefficients is provided. In Section III this 
expansion is used to find the critical temperature and specific heat of the Bose gas. In Section IV 
we summarize our main findings and explain how these can be used to extract information about 
extra dimensions present. 

II. PARTITION FUNCTION 

We consider a gas of ./V non-interacting bosons of mass M on a product manifold ^ = R 3 x ,yV 
of dimension D = d + 3, where the additional dimensions are modeled by a smooth, compact d- 
dimensional manifold JV with or without boundary djV . The dynamics of the quantum mechanical 
system we want to consider is described by the Schrodinger equation (here and in the following we 
set H = ks = 1) 

--i-A^ + F(x)^U n = £ n n , (2.1) 



2M 

where the coefficients E n represent the energy levels and we choose V(x) to be a 3-dimensional 
anisotropic harmonic oscillator potential 

V(x) = — (lu iX 2 + aj 2 y 2 + u 3 z 2 ) , (2.2) 
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which is the relevant choice for recent experiments. As is well known, equation (|2.ip can be solved 
by separation of variables, and the spectrum is found to be 

E a ,i = Xi + ^Uk (n k + -j , (2.3) 
k=i ^ ' 

with rik E No and where Aj denotes the eigenvalues of the Laplace operator A iy y on the manifold 

jV , so 

- Ajspi = Xnpi . (2.4) 

In order to study the thermodynamical properties of this system, the relevant object is the 
partition function or grand canonical potential 

q=-J2J2 ln [ 1 - e ~ 0iEn ' i ~^) ' ( 2 - 5 ) 
i n 

where fi denotes the chemical potential and we have introduced the standard notation (3 = 1/T. 
Since the ground state is of particular importance in the analysis of Bose-Einstein condensation, 
we separate its contribution from the rest of the series in (|2.5p and we expand the logarithm to 
obtain 

m=l i n 

where the prime indicates the omission of the ground state contribution and, denoting by go the 
degeneracy of the lowest eigenvalue Eq, 



Qo = -go In 



1 - e 



-p(E -n) 



(2.7) 



Mellin-Barnes integral representation [7J, 



The exponential that appears in the partition function (J2T6J) can be dealt with by exploiting a 

3,Q 



to obtain 



m=l _ i n 

r— too 

with r chosen in such a way that all the poles of the integrand lie to the left of the contour, and 
we have introduced the critical chemical potential [i c = Eq, since we are considering an ideal gas 
of bosons. 

After these manipulations we notice that the spectral ("-function 

cc«) = E'E' ^ - ( 2 - 9 ) 
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makes an appearance. The remaining sum over m in (|2.8j) can be written in terms of polylogarithmic 
functions [19j ] 



OO 7 
X 



1=1 



such that we are able to cast the partition function (|2.8p in the form 27h30I| 



r+ioo 



Q = q ° ~^Ti J d « r («)/ 3 ~ QL Wi(e- /3( " ) )c(a). (2.11) 

i — ioo 

The integral representation (|2,lip is particularly suitable for an asymptotic expansion as (3 — > 
0. Although for the phenomenon of Bose-Einstein condensation the low temperature regime is 
the pertinent one, it is clear for example from the references [23I . that the high-temperature 
expansion remains valid up to the condensation point and it can be used to determine the critical 
temperature at which condensation occurs. 

The high-temperature expansion of (|2.1ip is obtained by shifting the contour to the left. In 
that process we pick up contributions from poles of the spectral ^-function C( a )- The integral can 
then be computed by applying the Cauchy residue theorem. The position of the poles and the 
corresponding residues are found by exploiting the intimate connection of C( a ) with its heat kernel 

K< ^> = E'E^ - *^ - ^ ■ ( 2 - 12 ) 

i n 

The heat kernel in (I2.12D can be factorized as K(t) = Kii(t)Kj^(t), where Kn(t) is the heat 
kernel associated with the spectrum of the anisotropic harmonic oscillator and K,j/(t) is the heat 
kernel for the Laplacian on the manifold JV modified by the critical chemical potential which 
acts as a constant negative potential in (j2. 12j) . The small-t asymptotic expansion of (|2. 121) which 
encodes the residues of C( a ) can De obtained from the one for Kj{(t) and Kjf(t). The small- 
t expansion for the harmonic oscillator part is trivially obtained because the heat-kernel can be 
computed in closed form by observing that the sums are simple infinite geometric series. For K.y (t) , 



18, 



34 



under the made assumptions on ,yV, the small-t expansion is well known and reads 

^ 00 

y ' j=0,l/2,l,... 

with the heat-kernel coefficients s/^ of the modified Laplacian on jY . 
Combining the expansions of each factor, in detail we have 

, 00 /[fc/2] \ 

K(t) = L— -r ^ E GWif-i ^ > ( 2 - 13 ) 

wiw 2 w 3 (47r)2 k=Q y i=Q 2 J 
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where [x] represents the integer part of x, are the heat kernel coefficients on the manifold jV 
as given above, and the C[(uj), which come from the small-i expansion of K}j(t), have the form 



l n 

cm = (-iy £ £ f^j^w , 

j\{n- -n)\ 



(2.14) 



n=0 j=pO- 

with i?j denoting the Bernoulli numbers [19|]. From the knowledge of the asymptotic expansion 
of the heat kernel (|2. 13j) . one can show that the rightmost poles of Q{a) are located at = 
[d — k)/2 + 3, with k = 0, . . . , (d + 5) [18]. Their residues are given by 

[fc/2] 



ResC(« fc ) = T ^ CiM^f-i 



(2.15) 



Taking into account the first three rightmost poles of C( a )i we obtain the following asymptotic 
expansion for the partition function valid for j3 — > 



1 d+6 \ 



-/3(Mc-M) 



— ^ h/3 ^ 2 



-/3(p c -fJ.) 



1/2 



(2.16) 



+ 



fd+4\ 

/?-(— JLi 



rf+6 
2 



1 



+ (wi + o; 3 ) ) +0 (/?" 



1 



d+3 



The particle number N is of specific importance in the analysis of Bose-Einstein condensation since 
it is used in order to define and compute the critical temperature. It is well known that ./V can be 
obtained from q as 

1 dq 



N 



(2.17) 



with the derivative evaluated at fixed temperature and volume. The last remark, together with 
the result (|2.16p . allows us to write 



d+J>N 



N = N + /TV— JLid+6 



,-/3(mc-m) 







( d+5 \ 

+ P y 2 JLid+s 

WIW2W3 2 



3 -/5(Mc-M) 



'1/2 



W1W2W3 



(2.18) 



d + 4 



-/8(a*o-m) 



— ^- f ^ + + ^2 + w 3 ) ) + o ( r 



d+3 
2 



The partition function q also provides the energy of the system through the relation 

U = 



d [id 
"dp + ~pdji 



and it reads 



U = U + 



d + 6 



I d+8 \ 



d + S 
2 



3 -/3(Mc-M) 



W1W2W3 



+ 



d + 5 



f d+7 \ 

/3 ^ 2 JLid+7 



(2.19) 



+ 



d + 4 



( d+6 \ 
p v 2 /Lid+6 



d-l. 
2 



/ d+6\ 

+ /i c/ g-CT-JLi 

/ d+4 \ 

+ /x c /3 v 2 ) Li d +4 



-P(jic-lt) 







" ( *ff + o^(^l + ^2 + W3) 

9 l/2 



( d+5 \ 

+ M c /3 ^2 /Lid+5 
W1W2W3 2 



^1(^2^3 



1 



W1W2W3 



1 



^ + ^cf^i + ^2 + w 3 ) ) + O ( /3 



d+3 



(2.20) 
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This expression for the energy will be used, in what follows, in order to compute the specific heat 
of the Bose gas. 



III. CRITICAL TEMPERATURE AND SPECIFIC HEAT 

The critical temperature T c = 1//3 C , at which the condensate starts to appear, is obtained 
from the particle number N in (12.18|) by setting No = 0. When the temperature of the system 
approaches the critical temperature we have \i ~ fj, c and, hence, we can use the Taylor expansion, 
valid for n > 2, 



Li n (0 = ( R (n) - x( R (n - 1) + 0(x 2 ) . 



(3.1) 



In these circumstances, the critical temperature is approximately defined according to the relation 



n = /t-Wch ( 

and it can be found to be 



d + 6 
2 

d + 4 



1 



d + 5 



1/2 



*i + (wi + ^2 + oj 3 ) ) + O [p- 



d+3 



(3.2) 



to 



1 



2c« m ^ 



/2 



(d + 6)Ko; 2 a;3)^Ci?( i T) 



d+6\ d+6 



) d+6 



(3.3) 



iidJT^O and 



x _ 2 C* (jg) + + ^2 + o; 3 )] ^ 



when djV = 0. We would like to point out that in fj3.3j) and (|3.4p we have defined, for brevity, 



d+4 
<Jf\ d+6 



(3.4) 



2V 



Cr (ht; .-'0 



2 

d+6 



(3.5) 



Equation (|3.4|) shows in detail how the critical temperature correction due to the finite number 
of particles depends on properties of the extra dimensions and, of course, the frequencies of the 
rarmonic oscillator potential. In some detail, the leading heat kernel coefficients are well known 
18] and we have that equals the volume of JV , <^{/2 i s proportional to the volume of the 
boundary of Jf , and higher order terms involve curvature tensors of the manifold JV and its 
boundary. Thus a measurement of the critical temperature as a function of the finite particle 
number N could reveal properties of the extra dimension like its number and size. 
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A further interesting thermodynamical quantity associated with the system is the specific heat. 
It is derived from the energy according to the relation 

C-%, (3.6, 

where the derivative is understood as performed by keeping both the number of particles and the 
volume constant. From (|2.20p and (j3.6|) . near the critical temperature T c , one has 

c = <"+ g ><" + 'V (y)o, ( d -¥) < + ' < I 7|(d+ V ^{. (^) < 
+ofrH, (3.7) 



when the manifold <yf satisfies the property d.jV ^ 0, and 

(d + 6)(d + 4) fi ±4 V fd+6\f .j, , 1 



+ - : 1 - p~^Kr ( ^ ) ( < + + w 2 + ^ 3 



4wia;2W3 \ 2 / \ 2 

^-rtV + ^ft^V^y+ofH-'h (3.8) 



45o(wiw 2 w 3 ) 2 V 2 

for cL^f = 0. The results that we have obtained for the critical temperature (I3.3p - (]3.4p and for the 
specific heat (|3.7p - (|3.8p are very general and hold for an arbitrary smooth, compact manifold . 
More explicit formulas can be obtained once the manifold jY has been specified. In this case only 
the knowledge of the first few heat kernel coefficients is necessary, which are well known for a wide 
variety of manifolds with and without boundary 



32, 



4|. 



IV. FINAL REMARKS 



In this work we have presented an analysis of Bose-Einstein condensation on product manifolds 
by utilizing ^-function regularization techniques. The method proves to be very efficient and close 
to the condensation temperature the general results for the critical temperature and the specific 
heat depend only on the first few readily available heat kernel coefficients associated with the base 
manifold Jf . It is important to notice that the dependence of the critical temperature on the finite 
number N of particles reveals geometrical and topological properties of the extra Kaluza-Klein 
dimensions; see equations (j3.3j) and ()3.4p . 

It would be of particular interest to study in more detail the cases in which the manifold JV 
is either a (f-dimensional sphere or a (i-dimensional torus. In these situations the spectrum \ is 
explicitly known and for the sphere one would be led to deal with Barnes ^-functions, while for the 



8 



torus Epstein ("-functions would be the relevant objects. We plan to investigate these cases in a 
future work. Furthermore, similar calculations should be done for string inspired models such that 
future experiments could possibly determine features of extra dimensions in these models too. 
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